Nanosecond molecular relaxations in lipid bilayers studied by high energy resolution 

neutron scattering and in-situ diffraction 
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We report a high energy-resolution neutron backscattering study to investigate slow motions on 
nanosecond time scales in highly oriented solid supported phospholipid bilayers of the model system 
DMPC -d54 (deuterated l,2-dimyristoyl-sn-glycero-3-phoshatidylcholine), hydrated with heavy wa- 
ter. Wave vector resolved quasi-elastic neutron scattering (QENS) is used to determine relaxation 
times t, which can be associated with different molecular components, i.e., the lipid acyl chains and 
the interstitial water molecules in the different phases of the model membrane system. The inelastic 
data are complemented both by energy resolved and energy integrated in-situ diffraction. From a 
combined analysis of the inelastic data in the energy and time domain, the respective character of 
the relaxation, i.e., the exponent of the exponential decay is also determined. From this analysis 
we quantify two relaxation processes. We associate the fast relaxation with translational diffusion 
of lipid and water molecules while the slow process likely stems from collective dynamics. 

PACS numbers: 87.14.Cc, 87.16.Dg, 83.85.Hf, 83.10.Mj 



INTRODUCTION 

Inelastic scattering techniques give a wave vector re- 
solved access to dynamical properties, i.e., excitation fre- 
quencies or relaxation rates can be quantified at different 
internal length scales. This is important to associate the 
dynamics with individual modes of motion. A unique 
advantage lies in the simultaneous access to both struc- 
tural and dynamical properties by measuring the scat- 
tering S(Q,cj) as a function of wave vector transfer Q 
and energy transfer u. We hereby report on /LteV energy 
resolved spectra in the model membrane system DMPC- 
d54, achieved by the neutron backscattering technique. 
By analyzing the respective Q dependence, we simulta- 
neously probe low energetic density fluctuations of lipid 
acyl chains and interstitial water molecules, i.e., the wa- 
ter layer in between the stacked membranes. 

The spectrum of fluctuations in biomimetic and biolog- 
ical membranes covers a large range of time and length 
scales HEHSSHHIiH, ranging from the long 
wavelength undulation and bending modes of the bilaycr 
with typical relaxation times of nanoseconds and lateral 
length scales of several hundred lipid molecules to the 
short wavelength density fluctuations in the picosecond 
range on nearest neighbor distances of lipid molecules. 
Local dynamics in lipid bilayers, i.e., dynamics of in- 
dividual lipid molecules as vibration, rotation, libration 
(hindered rotation) and diffusion, has been investigated 
by, e.g., incoherent neutron scattering Q, 0,1110, HI and 
nuclear magnetic resonance ^| to determine the 
short wavelength translational and rotational diffusion 
constant. Collective undulation modes have been investi- 
gated using neutron spin-echo spectrometers 
and dynamical light scattering [lj,lla,lla| ■ Only recently, 
the first inelastic scattering experiments in phospholipid 



bilayers to determine collective motions of the lipid acyl 
chains and in particular the short wavelength dispersion 
relation have been performed using inelastic x-ray [T^ | 
and neutron 01 scattering techniques. While here fast 
propagating sound modes in the picosecond time range 
have been quantified, the present paper deals with slow 
nanosecond relaxation times on length scales of nearest 
neighbor distances of phospholipid acyl chains and wa- 
ter molecules, i.e., the slow dynamics of melting (diffu- 
sion) and collective movements (most likely undulations) 
of the lipid and water backbone. We have selected the 
neutron backscattering technique for this study since the 
dynamical modes at high q\\ (the in plane component of 
the scattering vector Q) are too fast to be accessed by 
x-ray photon correlation spectroscopy (XPCS) and the 
lateral length scales are too small to be resolved by dy- 
namic light scattering (DLS) or neutron spin-echo tech- 
nique (NSE). Note that the feasibility of wave vector re- 
solved backscattering experiments in oriented membrane 
systems to study freezing of lipid acyl chains and wa- 
ter molecules has been established only recently [l9j . 
In this work we present the unprecedented determina- 
tion of relaxation rates from Q-resolved backscattering 
QENS data in phospholipid bilayers, combined with in- 
situ diffraction to unambiguously assign the relaxation 
times to the respective phases of the model membranes. 
These measurements are carried out in aligned phases to 
preserve the unique identification of modes on the basis 
of the parallel and perpendicular components q\ \ , the lat- 
eral momentum transfer in the plane of the bilayers, and 
q z , the reflectivity, of the scattering vector Q. 
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EXPERIMENTAL 

In the case of single membranes the inelastic neu- 
tron scattering signal is by far not sufficient for a 
quantitative study of the inelastic scattering. Partially 
(acyl chain) deuterated DMPC-d54 (1,2-dimyristoyl-sn- 
glycero-3-phoshatidylcholinc) was obtained from Avanti 
Polar Lipids. Highly oriented multi lamellar membrane 
stacks of several thousands of lipid bilayers were prepared 
by spreading lipid solution of typically 25mg/ml lipid in 
triffiioroethylene/chloroform (1:1) on 2" silicon wafers, 
followed by subse que nt drying in vacuum and hydration 
from D2O vapor |2fj, resulting in a structure of smec- 
tic A symmetry. Twenty such wafers separated by small 
air gaps were combined and aligned with respect to each 
other to create a " sandwich sample" consisting of several 
thousands of highly oriented lipid bilayers (total mosaic- 
ity about 0.5°), with a total mass of about 400 mg of 
deuterated DMPC. The sample was mounted in a her- 
metically sealed aluminium container within a cryostat 
and hydrated from D2O vapor. Saturation of the vapor 
in the voids around the lipids was assured by placing a 
piece of pulp soaked in D2O within the sealed sample con- 
tainer. The pulp was shielded by Cadmium to exclude 
any parasitic contribution to the scattering. The hydra- 
tion was not controlled but we allowed the sample to 
equilibrate for 10 h at room temperature before the mea- 
surements. The large beam divergence as compared to 
reflectometers did not allow to determine the correspond- 
ing d z spacing simultaneously with sufficient accuracy to, 
e.g., determine the swelling state of the membrane stack. 
We note that the absolute level of hydration is therefore 
not precisely known but the temperature of the main 
transition agrees quite well with literature values. 

The experiment was carried out at the cold neutron 
backscattering spectrometer IN16 jilj at the Institut 
Laue-Langevin (ILL) in its standard setup with Si(lll) 
monochromator and analyzer crystals corresponding to 
an incident and analyzed neutron energy of 2.08 meV 



FIG. 1: (Color online), (a) Energy integrated diffraction for 
temperatures 250 K<T<310 K with a temperature resolution 
of AT «0.3 K (normalized to acquisition time). q\\ is the in- 
plane component of the scattering vector Q. The phase bound- 
aries for gel (L/3), ripple (P/3') and fluid phase (La) as defined 
by structural changes, i.e., the change of peak position and 
width are marked by the solid white lines, (b) 'True' elastic 
scattering, measured with an energy resolution of 0.9 /j,eV. The 
phase transitions are also visible in the energy resolved diffrac- 
tion data by following the melting at the lipid acyl chain and 
water position. Data are normalized to monitor, detector ef- 
ficiency normalized to the scattering signal of a 2 mm thick 
Vanadium plate oriented at 135° with respect to the incom- 
ing beam. Please note that no absorption correction has been 
done. For geometrical reasons the diffraction detectors have a 
lower maximum qii than the backscattering detectors. 

(A=6.27 A) resulting in a high resolution in energy trans- 
fer of about 0.9 /j,eV FWHM. An energy transfer of 

— 15^eV< E < +15^eV can be scanned by varying the 
incident energy by Doppler-shifting the incident neutron 
energy through an adequate movement of the monochro- 
mator crystal. 20 detectors in exact backscattering with 
respect to the analyzer crystals covering an angular width 
of 6.5° each have been used. Note that the maximum of 
the static bulk water correlation peak at about q\\ ~2 A -1 
is not accessible because of simple geometrical reasons. 
The set-up is therefore more sensitive to detect amor- 
phous water with slightly larger nearest neighbor dis- 
tances, as it is expected to occur in the water layers of 
the stacked bilayers. A separate line of 160 energy in- 
tegrating diffraction detectors - normally mounted on a 
parallel circle inclined by 24.4° below the scattering plane 

- with an angular width of 1° each - allows to simulta- 
neously detect structural changes with a much higher q\ | 
resolution. 



RESULTS FROM NEUTRON DIFFRACTION 

We thus simultaneously performed two types of diffrac- 
tion measurements, namely energy-resolved and energy- 
integrated measurements. Note that the additional 
diffraction detectors very well integrate over low lying ex- 
citations and relaxations. Figure 2] (a) shows an energy 
integrated diffraction pattern for 1.1 A _1 < gii<1.85 A -1 . 
From the temperature dependence of the inter-acyl chain 
correlation peak at <7| | wl.4 A -1 we assign the transition 
temperatures of gel to ripple phase (L/3-P/3') to 285 K 
and the temperature of the main transition (P/?'-Lcv) to 
T=295 K. T/ u) =271 K marks the freezing temperature 
of the amorphous water layer in between the stacked 
membranes |l9j. The correlation peak gradually shifts 
to smaller q\\ values below Tf w until about 264 K. In the 
range 264 K<T<271 K water migrates out the bilaycr 
stacks to freeze as bulk water, as discussed in detail in 
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FIG. 2: (Color online). Inelastic scans at the lipid chain 
(gn ~1.4 A -1 ) and water (gn wl.9 A -1 ) position for tem- 
peratures T=250 K, 280 K and 310 K. Detectors have 
been grouped to increase counting statistics. The inter 
acyl chain correlation peak is measured by detectors 9- 
13, which cover a q\\ range of 1.20 A _1 < g N <1.60 A -1 . 
The water contribution is detected by detectors 17-20 
covering 1.78 A _1 < q||<1.94 A" 1 . As a model, up to 
two Lorentzian peak profiles (dashed lines) to describe 
the quasi elastic broadening and a Dirac function (dot- 
ted line) to describe the elastic intensity were assumed. 
This model was convoluted with the measured resolu- 
tion function obtained from a Vanadium standard. A 
flat background (dashed line) was subsequently added 
and the result was fitted (fit result: solid line) to the 
data (circles). A flat background may arise from fast 
processes far beyond the accessible energy window of the 
spectrometer. The thus obtained relaxation times are 
given in the figures. The calibration error of the energy 
scale is approx. 3%. 
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FIG. 3: (Color online). Corresponding Fourier trans- 
formed energy space data (power spectra) from Fig. |2] 
(points), divided by an analytical description of the 
Vanadium resolution. The thus obtained intermediate 
scattering function I(qn,t)/I(qi\, 0) is shown together 
with the fits (blue solid line) of two exponential decays 

/(g||,t)/7(9ll,0) = {Ao-A^eM- (*Ao( 9 | |))^°] + yi + 

(Ai — yi)exp[— (£/i"i(q||))' 31 ]. The values for to.i have 
been fixed from the energy scans and only the exponents 
/3o,i and the amplitudes Ao,i have been fitted. The thus 
obtained values for /3o are given in the figures. The res- 
olution is also plotted (green line). While the data at 
250 K are well described by one exponential decay, two 
decays are needed at higher temperatures. 



Ref. 0|- The energy resolved pattern in Fig.Q](b) sheds 
light on the melting process of lipid acyl chains and water 
and perfectly reproduces the results previously reported 
in [19( (the difference being the larger number of detec- 
tors). Only structures which are static with respect to 
the energy resolution of 0.9 /xeV, which correspond to 
motions slower than approx. 4 ns, are detected as elastic 
when using the energy resolving backscattering detectors. 
Changes in this intensity may arise from either structural 
or dynamical changes, i.e., shifts of correlation peaks or 
freezing or melting of dynamical modes. Elastic intensity 
at the water detectors rises at T<271 K; at the acyl chain 
position there is an elastic contribution below 295 K, the 
temperature of the main transition. A further increase 
occurs when entering the gel (L/3) phase indicating the 
high degree of dynamics in the ripple (P/3') phase. Note 
that there is very little elastic contribution in the fluid 
phases pointing out that the intensity usually measured 



in diffraction experiments stems predominantly from low 
lying fluid excitations and quasi elastic scattering. 

INELASTIC SCATTERING 

Based on the diffraction patterns we have recorded in- 
clastic scans at temperatures of T=250 K, 280 K, 292 K 
and at 310 K, far in the fluid La phase of the model 
membranes. Additional scans in the regime of critical 
swelling (T=296, 297, 299 and 303 K) have also been 
taken. The results for the acyl chain (qu wl.4 A -1 ) 
and the water position (qii ml.9 A -1 ) are depicted in 
Fig. At 250 K, both g| positions are well fitted by 
a single Lorentzian, i.e., a single relaxation process with 
an energy width (FWHM) Aw and corresponding relax- 
ation time t = 2n/Aui. For temperatures above 280 K 
for the water and above 292 K for the acyl chains two 
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FIG. 4: (Color online), (a) Relaxation times at the lipid 
chain and the water position for all measured temperatures 
as determined from fits to the quasi elastic data in Fig. [5] (b) 
Exponents of the exponential decay as determined from fits of 
the intermediate scattering function from Fig. [3] Solid lines 
are guides to the eye. 



Lorentzians are needed to describe the data. Although 
the Lorentzians describe the energy spectra, small de- 
viations from the fits are visible pointing to deviations 
from single exponential relaxations. A detailed analysis 
in energy space is difficult and would mean to add more 
Lorentzian or mix Lorentzian and Gaussian peak profiles 
with the risk of being physically meaningless. To fur- 
ther characterize the relaxation processes the data have 
been Fourier transformed (using the Matlab FFT algo- 
rithm) to determine the exponent (3 of the corresponding 
exponential decays in time domain, exp [— (t/r) 13 }, from 
fits to the intermediate scattering function I(qu,t), as 
shown in Fig. [3] With increasing temperature there is 
a relaxation step moving into the time window at the 
two <gr 1 1 positions, respectively. Because of the limited dy- 
namical range of 15 fJ,eV, the shortest time accessible is 
about 0.2 ns and the faster process is already partially 
out of the accessible time window, i.e., the quasielastic 
broadening is slightly broader than the energy window. 
No exponent could therefore be reliably determined for 
the fast process. The resulting relaxation times to,i and 
the corresponding exponent of the slow process, /?o, are 
depicted in Fig. ^ (a) and (b) for all measured temper- 
atures. To gain reliable fitting parameters for r and /3, 
the relaxation times ro,i have been determined from the 
energy spectra S(q\\,u) in Fig. |2] because this gives the 
highest accuracy and stability. The r values have then 
been fixed for the fits of the exponential decay in the time 
domain. f3 is determined from the fits to I{q\\,t) in Fig. 01 
The two processes at the two q\ \ positions show distinctly 
different relaxation times and temperature dependencies. 



While the fast process has constant times of about 0.1 ns, 
the slow process starts at about 1 ns at 250 K and be- 
comes faster as approaching the main phase transition at 
295 K. Note that because of the limited dynamical range 
the accuracy in determining the fast relaxation time T\ 
is also limited and we can not exclude a temperature de- 
pendence of the fast process. While for the fast branches 
no exponent can be determined for the same reason, the 
slow processes show a stretched exponential character at 
low temperatures and turn to single exponential when 
the corresponding correlations melt at the lipid (main 
transition) and water melting temperature, respectively 



DISCUSSION 

From the inelastic neutron data we quantify two relax- 
ation times. We argue that the corresponding processes 
stem from collective motions of the lipid acyl chains 
and interstitial water molecules, respectively, rather than 
from local dynamics. By selective deuteration of the 
chains and hydration from D2O vapor, coherent motions 
are strongly enhanced over other contributions to the in- 
elastic scattering cross section The fast process T\ 
disappears below the freezing temperature of the lipids 
(P/3' phase) respective water molecules (T<T/ W =271 K) 
and might therefore be due to (collective) translational 
diffusion of lipids and water. The slow process with re- 
laxation times To is measured down to the lowest temper- 
atures of T=250 K and most likely probes collective dy- 
namics, i.e., neighboring lipid and water molecules partic- 
ipating in slow mesoscopic dynamics. However, to further 
characterize this process, q\\ dependent data are needed 
to determine the corresponding dispersion relation and 
compare to theoretical models. It is speculated in the 
literature that the minimum in the dispersion relation 
of the collective fast propagating modes at the nearest 
neighbor distances of lipid acyl chains (as determined in 
Rcfs. H3 and [3) is related to transport phenomena 
within and across the bilayers. The relaxation found in 
the present paper on the same length but distinctly dif- 
ferent time scale might be relevant to establish such a 
model of phonon assisted diffusion in membranes, which 
would be of particular interest in membrane biophysics 
and biotechnology applications. 

The slow relaxations at the acyl chain and the wa- 
ter position both show a decrease of relaxation times 
towards the main transition of the phospholipids. The 
relaxation time of the fast branch with r «0.1 ns that 
we attribute to translational diffusion is temperature in- 
dependent within the experimental accuracy, which is 
slightly limited by the maximum achievable energy trans- 
fer. Note that the diffusion time measured here can be in- 
terpreted as time for hopping processes to nearest neigh- 
bor sites. Such a process needs a free neighboring site but 
does not necessarily give information about the mobility 
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of these holes (free surface or volume) what would be 
important for macroscopic diffusion. The corresponding 
time scales and diffusion constants might therefore be 
distinctly different. From measurements with a higher 
dynamical range, the energy barrier of this excitation 
could probably be determined from the temperature de- 
pendent hopping times using Arrhcnius or Vogel-Fulcher 
laws. 

The (3 exponents in Fig. 01 (b) start at values (3 < 1 at 
low temperatures. Structural inhomogeneitics and het- 
erogeneous interactions obviously lead to a local relax- 
ation dynamics and to stretched (J3 < 1) exponentials. 
When entering the fluid phase at T m , the lipid relax- 
ations (q|| w 1.4 A" 1 ) turn into single exponential, as 
can be expected for a diffusive, fluid like motion of the 
particles. At the water position (qn rs 1.9 A -1 ), the 
transition from stretched to single exponential occurs be- 
tween 250 K and 280 K, most likely at the water freezing 
or melting temperature T/ u ,=271 K, as speculated and 
indicated in the Figure. A striking feature is that in the 
range of critical swelling, the water relaxations become 
again stretched what might be related to the expansion 
of the water layer leading to the well known anomalous 
swelling of phospholipid bilayers close to the main tran- 
sition [22j. 

CONCLUSION 

In conclusion, using the neutron backscattering tech- 
nique, we have determined q\\ -dependent relaxation times 
at different temperatures in the model membrane system 
DMPC that we attribute to lipid acyl chains and intersti- 
tial water molecules. In-situ diffraction allowed to clearly 
assign the temperatures to the different phases of the bi- 
layers. A combined data analysis in the energy and time 
domain allowed to quantify relaxation times of two pro- 
cesses and determine the exponents (3 of the exponential 
decays. The fast process is attributed to diffusion while 
the slow process likely stems from collective dynamics. 
We find stretched exponential relaxations at low tem- 
peratures. The relaxations turn into single exponential 
above the corresponding melting temperatures of lipid 
and water molecules. 

Future experiments, which will include a selective 
deuteration of the acyl-chains and the membrane water, 
respectively, will allow to mask different types of mobility 
and hopefully deduce complete dispersion relations of the 
different molecular components in the different phases of 
the phospholipid bilayers. Compared to NSE and XPCS, 
large Q-values can easily be obtained. The current draw- 
back of the technique is the limited dynamical range on 
high-flux spectrometers which hopefully will be overcome 
in next generation neutron sources and backscattering 
spectrometers. Our study also points to the wealth of 
new information that can be explored with the membrane 



backscattering technique. 
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